Correlated tunneling in intramolecular carbon nanotube quantum dots 
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We investigate correlated electronic transport in single-walled carbon nanotubes with two in- 
tramolecular tunneling barriers. We suggest that below a characteristic temperature the long range 
nature of the Coulomb interaction becomes crucial to determine the temperature dependence of the 
maximum G max of the conductance peak. Correlated sequential tunneling dominates transport 
yielding the power-law G max oc T a °' ld ~<" ld ~ 1 , typical for tunneling between the ends of two Luttinger 
liquids. Our predictions are in agreement with recent measurements. 
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Electronic correlations have been predicted to domi- 
nate the characteristic features in single-walled metal- 
lic carbon nanotubes (SWNTs) jl], 1, which has re- 
cently been observed in experiments [|, ||, |). The one- 
dimensional nature of the electronic conduction bands 
reveals itself in typical Luttinger liquid rather than 
Fermi liquid behavior. The manipulation of individual 
nanotubes with an atomic force microscope permits the 
creation of intra-tube buckles acting as tunneling barri- 
ers || . Recently, SWNTs with two intramolecular buck- 
les have been reported to behave as room-temperature 
single-electron transistors 0. At low temperatures, the 
thermal electronic energies are smaller than the level sep- 
aration between the discrete energy states of the island 
and tunneling occurs via these discrete levels. Resonant 
tunneling in Luttinger liquids has been investigated the- 
oretically by many authors f| g, fuj [Q|> @. The one- 
dimensional nature of the correlated electrons is responsi- 
ble for the differences to the quantum Coulomb blockade 
theory for conventional, e.g., semiconducting quantum 
dots JT3| . Varying the gate voltage results in a sequence 
of conductance peaks. In the (uncorrelated) sequential 
tunneling (UST) approximation the temperature depen- 
dence of the maxima of those peaks follows the power-law 
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with 



being the density of states exponent for tun- 



neling into the end of a Luttinger liquid. However, recent 
experiments M suggest a different power-law, 
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In this Letter, we propose that a novel tunneling mech- 
anism, correlated sequential tunneling (CST), gives rise 
to the power law (g). It originates from the finite range 
nature of the Coulomb interaction in SWNTs and re- 
places conventional uncorrelated sequential tunneling. It 
dominates resonant transport at low temperatures and 
strong interactions. CST leads to a renormalization of 



the intrinsic linewidth of the resonance, which is reflected 
in an increase of G max with increasing temperature. In 
contrast, UST would predict the opposite behavior Jn4| . 
A good agreement with the experimental results is 
found. 

We describe an individual metallic SWNT with two 
buckles by the Hamiltonian H = Ho + Hb + H cyLt . Here, 
Hq characterizes the one-dimensional homogeneous wire 
including the finite-range electronic interaction. Metallic 
carbon nanotubes possess two gap-less one-dimensional 
bands j|, |) with Fermi velocity vf — 8 x 10 5 m/s §, 
which dominate the low-energy physics. In the bosonized 
representation, one finds 
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+ - dx dx'd x d p (x)V{x - x')d x S p (x') , (3) 



containing three sectors of neutral excitations (a = 
a, Ap, Act) propagating with vf, and one sector of 
charged excitations {a = p). The charge sector is char- 
acterized by the nonlinear dispersion relation u)(k) = 
VF\k\[l + V(k)/TrhvF], where V(k) is the Fourier trans- 
form of the screened Coulomb potential projected onto 
the x direction V(x) |l7| . In our approach, the non- 
linear dispersion is crucial, since it modifies the low 
energy properties of the dynamics. The two tunnel- 
ing barriers at x = ±xq/2 are modeled by a scat- 
tering potential U(x) = CbE^^^ 1 + KX o/2) lead- 
ing to Hb = J dxll(x)p(x). Here, p(x) is the total 
electron density and backscattering has to be included. 
Upon introducing N a = 4= [§ a (?f) +tf a (-^)] and 

n a = [0„ (l) - Va (-f )] + ^S a , P , H B provides 
an eight-dimensional sinusoidal potential for the tunnel- 
ing dynamics associated with the variables N a and n a 
S. Physically, eN p /2 counts the unit charges trans- 



ferred through the dot, while 



is the total charge 



accumulated on the dot. Likewise, N a and n a are re- 
lated to the difference of spin fluctuations between both 
leads, and to the deviation from the mean value in 



2 



the dot, respectively. The analogous interpretation for 
Nao-, it-Act and NA P ,riAp holds. Finally, the gate volt- 
age Vg and the bias voltage V give rise to the term 
H cxt = -e[{V/2)Np + XV G n p ], where A = C G /C toU with 
Cg the gate capacitance and Ctot the total capacitance 
of gate and island. 

The relevant observable is the asymptotic tunneling 
current / = lim^oo ^(N p (t)), where (...) denotes the 
expectation value determined via the reduced density 
matrix (RDM) . The latter is obtained by tracing out ex- 
actly the Luttinger modes away from the barriers JL8]. 
The diagonal elements of the RDM are given by the con- 
ditional probabilities P(Nf , t; N;, to) of being in the final 
state Nf := ({iV a> f }, {n a ,i}) at time t, having started 
from the state Ni at time to- Then, we find 

I = \ lim V N Pti P(N f , t; N h t ) := e(T f - T b ) , (4) 



be evaluated analytically. In general, as in the zero-range 
case, V zr (x) = Vo5(x), the leads contribute an Ohmic 
spectrum at low frequency: lim^^o J+( w ) + J-( w ) := 
J£(cj) = u/g p , where g~ 2 = 1 + V(k = 0)/7rfru F . Theo- 
retical estimates as well as experiments show that typi- 
cal values for g p range between 0.18 and 0.26 in SWNTs 
[jl], 0, pi, 0], being well below g p = 1 for noninteracting 
electrons. 

As a consequence of the finite range interaction, a new 
energy scale arises which introduces two qualitatively 
new features: (i) The 5-shaped resonances for the zero- 
range limit turn into resonances of finite height and with 
an intrinsic line- width oc 1/gA, cf. Fig. |l|. (ii) The new 
parameter g^ 1 is defined from the low energy relation 



lim J+(w) — J-(uj) :— J A {u) = co/g A 



with r f / b being the total rates for forward/backward 
transfer of charge. In turn, after exploiting the detailed 
balance relation T b (y) = e~^ eV T f (V), the linear conduc- 
tance follows as G = limy^o e 2 /?r f . Here, (3 = l/(fcBT) 
is the inverse thermal energy. 

In general, the conditional probability P(Nf , t; Nj, to) 
can be evaluated exactly in terms of a multiple real 
time path integral involving forward N a (r) (n a (r)) and 
backward N' a {r') «(r')) paths The effect of 

the bath modes d a {x ^ ±xo/2,t) enters via the influ 
ence functionals J-+ = exp{^ Q $ + [7V a , N' a ]}, and T- = 
cx P{J2a ^~\ n a,n' a ] + $c[na,n' a }}, where the phases 



$ ± {q,q'} = f dr f dr'[q(r) - q' \r')} 



x{WUt ~ r')q(r') - W?(t - r')g'(r')}5) 

induce nonlocal-in-time correlations among different tun- 
neling transitions. The correlation functions read J18[ 



For ease of calculation, we assume a potential V(x) = 
Voa/2 exp(— a\x\) |l7]], with a being the inverse screening 
length (a — » oo yields the zero-range interaction limit 
VoS(x) and implies g^ 1 — ► 0). The numerical results for 
this potential are shown in Fig. [|. Note however that, 
as shown below, our main conclusions are independent of 
the precise shape of the potential. Finally, $c provides a 
quadratic local-in-time contribution proportional to the 
addition energies E a = AE a /2 for the neutral sectors, 
= and E p = K~(u> n = 0) for the charge plasmons. 

Of foremost importance is how to evaluate the for- 
ward rate T . In the large barrier limit, the dynamics 
is ruled by well separated tunneling events between the 
adjacent minima of the periodic potential in Hq. Each 
tunneling event induces the change N a — > N a ± 1, and 
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with the spectral densities 
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use 
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(AE a /H) J2m=i — ^±( m ))]/4 f° r the neutral sectors 
(a = a, Act, Ap). Here, we have introduced a cut-off 
frequency u) c , and Q + (m) = AE a (2m — l)/h, £1_ = 
AE a 2mjh. Finally, AE a = nhvp/xo is related to the 
energy quantization of the neutral plasmon modes in 
the dot. In 0, AE a > k^T also at room tempera- 
ture. For the charged sector the spectral density reads 



j p ± {uj) = 



Q + )/ttH, where 



[K ± K) = — 
nir 



~ dk l±CQ S (kx ) ^ (?) 



and w n are the Matsubara frequencies. For an arbitrary 
finite range potential the spectral densities J±(cj) cannot 



n a — > n a ± 1. In particular, N p — > N p ± 1 describes one 
discrete charge tunneling through one barrier from left to 
right (— ) or vice versa (+), while n p — > n p ± 1 describes 
a unit charge tunneling onto/out of the dot. Each tun- 
neling event through the left (right) barrier contributes 
a factor iAl(— ^Ar), where Al = Ar, is the associated 
tunneling amplitude. To proceed, we express the dou- 
ble path integrals over the paths N a ,N^ and n a ,n' a as 
single path integrals over paths along the states of the 
RDM in the (N a ,N' a )- and in the (n OJ nJJ-planes. Ac- 
cording to the boundary conditions, all paths start and 
end in diagonal states of the RDM, making intermittent 
visits to off-diagonal states, cf. Fig. ^. Due to the strong 
Coulomb interaction (g p < 1), the relevant paths are 
those for which the system is back to a diagonal state 
after two tunneling events pl|. We focus on the linear 
regime (V — > 0), so that either n p = n or n p = n + 1 
units of charge are allowed in the dot in the stationary 
state The current (||) is then expressed as a series 
expansion in the tunneling amplitudes A R / L . To lowest 
order, the forward rate is divergent. The divergent term 



3 



reads 



It reads at resonance (E K 



0) 



A 2 A 2 roo ( 3 \ 

rf.w = iim y ( JJdriJ e -^-^(-o-^(- 3 ) 

x cosh Af 3 cos ^f 3 [cLf(ri)cRf(r 3 ) + c R f (n)c L f (r 3 )], (9) 



where Tj denote the time intervals elapsed between the 
different tunneling events. Then, T2 is the time spent 
in the intermediate diagonal state. We use the notation 
c Ka (r) = cos(i?s(r) — E K:a r), where the energies are de- 
fined as E^{(n) — — fi(n + 1) — eV/2 and E R j(n + 1) = 
/i(n + 1) — eV/2. Here, the chemical potential reads + 
1) = £ p (n - n Q - C G V G /e + 1/2) + (3/2)£ CT (<r 2 +1 - a 2 n ), 
with no = Ak-pXo/Tt being the mean electron number on 
the dot. Moreover, a n = (±1) for an even (odd) num- 
ber of spins in the dot. The correlations are encapsu- 
lated in 5s/ a (t) and -Rs/a( t ) being the real and imag- 
inary parts, respectively, of VFs/A = Y^ a 0^+ ^ W-)- 
The long range nature of the Coulomb interaction is re- 
flected in the correlators Afg, and TZf^, inducing dipole- 
dipole correlations between the tunneling events across 
the left and right barrier, cf. Fig. |. To be definite, Af 3 = 
S a (Ti+T 2 +T 3 ) + S a (T 2 ) -5 q (ti+T2)-5 q (t2+t 3 ), and 

K% = R a {Tl+T 2 +T ? ,)-R a {T 2 +T 3 ,)+R a {T l +T2)-R a {T2)- 

Here, a = A if the interaction is between tunneling 
through different barriers, while a = £ if events through 
the same barrier are involved. Due to (^), the correla- 
tors associated to behave as 5 a (t) — > itk^TT / g&h, 
const, at long times. Hence, A^g P — * and 
The functions 5 A and R% 
associated to the neutral sectors are purely periodic with 
frequency AE a /h. 

To cure the divergence of the fourth order rate, we 
have to sum up divergent irreducible terms of higher or- 
der. Therefore, appropriate approximation schemes arc 
required. We make the following approximations: (i) We 
only consider higher order paths which yield the UST re- 
sult for vanishing interaction range, (ii) We neglect the 
dipole-dipole correlations A?- due to g p <C 1, but keep 
the correlations A A among the inner dipoles up to linear 
order, hi) We consider a large cut-off to c (scaling limit), 
and neglect the correlations Ttfj ,f '• The sum over the 
so obtained higher order irreducible paths yields for the 
total rate T f = ^2^L 2 r f, ( 2 ") a finite result. It can be ob- 



2 a 2 r°° 



J\dT t 



and R P A 
1Zi 3 p — > when t 2 — > oo. 



tained from (B) upon replacing Afg — > 3A 13 a , 7Z 13 P = 0, 
and A — » A + 7. Here 7 can be interpreted as a linewidth, 
i.e., 7 = T tot + 7 A , where T tot = T Rf + T u + T Rb + T Lb , 
and T Ka = (A 2 J2) f Q °° dTe~ s ^c Ka {r) are the incoher- 
ent rates for forward/backward tunneling through the 
right/left barrier (E L j(n) = -E L , h (n + 1), E Rth (n) = 
—E Ry {(n+l)). In contrast, j A is the lowest order correc- 
tion to the linewidth due to the finite range interaction. 



e -S s (r 1 )-5' s (r3) e -r t °V 2 



x siniJa(n) sin^s^Af/ coshSAf/ . (10) 

The incoherent rates yield at resonance to the power- 
law T K ' a oc T" ond , where a cn d = l/ffcff — 1, and 



l/ 5cff = (l/g p + 3)/4. In the limit A^ p = 
one recovers the known UST rate jy], |l2| T l = 
r Rf r Lf/ r tot K T a cnd; which imp i ies the p OW e r l aw 

(|l|) for G max . While r tot decreases with decreasing 
temperature, 7 A remains constant. Hence, depending 
on g e g a cross-over temperature T* can be identified, 
with k B T*/huj c = (7 A /d goff ) 1 ~ 1 / 9rff /27r with d ga{{ = 
A 2 |r(l/2 5cff )| 2 /( Wc r(l/ 5cff ) C ). For T < T*, the power- 
law (||) follows. Due to the nonlinearity of the spectral 
densities J£,^(uj), we are able to compute the triple in- 
tegral numerically only for rather large inverse screening 
lengths a, since then the correlators W^ a (t) are well 
fitted by those of a damped harmonic oscillator. The 
result for a — 200xq , yielding g^ 1 — 0.087, is shown 
in the right inset of Fig. 1 for g p — 0.23. Turning to 
the experiment in m we should a tyical screening length 
a -1 « 3 — 5nm. For a dot length xo — 25nm it corre- 
sponds to values a ~ 10xq . However, due to the strong 
asymmetry of the peaks of the spectral densities J^/a 
at a ~ 10xq , a numerical evaluation of 7 A is no longer 
possible. In the following, we assume a = 200:Eq 1 . 

We now compare the outcomes of (|^) with the experi- 
mental results Q . Induced buckles on a SWNT define an 
island of ~ 25 nm length. We include the contact influ- 
ence via the Landauer formula, in which the total tube 
conductance is convoluted with the Fermi distribution 
from the Au-leads. With typical Luttinger parameters in 
the range g p w 0.18 — 0.26, the agreement between the 
measured data and the theoretical prediction is very good 
upon choosing g p = 0.23, leading to g e g — 0.54. The 
comparison is depicted in Fig. ^ together with the power 
law (||). Note that cotunneling events, which should dom- 
inate the tails of a conductance peak [|[ |l^, are not 
considered here. The absolute value of G max is fitted 
by the parameter Al/r, (we have chosen fao c = 9AE a ). 
Importantly, we find a very good agreement between the- 
oretical and experimental results for k^T -c Ai? CT with 
the only one relevant free parameter g p . 

In the regime k^T ~ A_E CT , we observe a deviation from 
the power law due to the contribution of excited plasmon 
levels which become thermally accessible, as well as due 
to a crossover to the UST limit. The result of our simple 
model with two 5— shaped tunneling barriers are fitted 
to the data. We find x = 71 nm (AE a = 23 meV) 
which agrees with the measured value of xo — 25 nm. A 
smaller value for xo only shifts the onset of the deviation 
from the power-law to higher temperatures. Indeed, the 
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measured height profile of the buckle region suggests that 
the real scattering potential is more complicated. How- 
ever, our central result, namely the algebraic behavior 
at low temperature, is not affected by the precise shape 
of the impurity potentials. We also fit the theoretically 
predicted line-shape of a conductance peak to the mea- 
sured data as displayed in the left inset of Fig. |3| and 
find A = 0.25, which coincides with the experimentally 
determined value. We measure the width of a calculated 
conductance peak and find that it increases linearly with 
increasing temperature, in agreement with the measured 
data (right inset of Fig. |J). To underpin the algebraic 
behavior, we also show in Fig. || the conductance G* ob- 
tained by integrating the linear conductance over the gate 
voltage. Due to the linear increase of the peak width with 
temperature, the exponent of the power law is increased 
by one, in clear agreement with the experiment. 

In summary, we have shown that the finite range of 
Coulomb interaction can play a crucial role in the reso- 
nant tunneling in one-dimensional quantum dots. 
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FIG. 1: Spectral densities J^ a (oj)/lj for finite interaction 
range (a — 200a; ( J' 1 , g p = 0.23) with uj a = AE a /h. Left in- 
set: Enlargement of the low frequency part. The important 
parameter is 1/<?a = lim^^o J /lu. Right inset: Contri- 
butions to the effective line- width r tot and 7 A , with the latter 
being originated from the finite range of the interaction (see 
text). 
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FIG. 3: Conductance maximum G max (right) and integrated 
conductance G* (left) as a function of temperature. Solid line: 
theoretical prediction with A_E CT = 23 meV and interaction 
strength g p = 0.23.B: experimental results [7], dashed lines: 
power laws, see Eq. (||). Left inset: A single conductance peak 
for varying gate voltage for T =4.4 K, 20 K, 30 K. Right inset: 
Peak width versus temperature. 



